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Abstract. We calculate the redshift-space power spectrum of the Sloan Digital Sky Survey (SDSS) Luminous Red Galaxy 
(LRG) sample, finding evidence for a full series of acoustic features down to the scales of ~ 0.25 /jMpc"'. This corresponds up 
to the 9. peak in the CMB angular power spectrum. The acoustic scale derived, (107.8 ± 4.3) A"' Mpc (k < 0.2 h Mpc"'), agrees 
very well with the "concordance" model prediction and also with the one determined via the analysis of the spatial two-point 
correlation function by Eisenstein et al. (2005). This is not only an independent confirmation of Eisenstein et al. (2005) results 
made with different methods and software but also, according to our knowledge, the first determination of the power spectrum 
of the SDSS LRG sample. By calculating the two-point correlation function using the smooth cubic spline model fitted to the 
observed bandpowers and comparing with the results of the direct determination we demonstrate the consistency of our results. 
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1. Introduction 

In the beginning of 1970's it was already realized that acoustic 
waves in the tightly coupled baryon-photon fluid prior to the 
epoch of recombination will lead to the characteristic maxima 
and minima in the post-recombination matter power spectrum. 
The same mechanism is also responsible for the prominent 
peak structure in the CMB angular power spectrum (Sunyaev 
& Zeldovich 1970; Peebles & Yu 1970; Doroshkevich et al. 
1978). The scale of these features reflects the size of the sound 
horizon, which itself is fully determined given the physical den- 
sities Q.t,h^ and Q,„/i^. The acoustic horizon can be calibrated 
using the CMB data, thus turning it into a standard ruler which 
can be used to carry out standard cosmological tests. For exam- 
ple, if we are able to measure the redshift and angular intervals 
corresponding to the physically known acoustic scale in the 
matter power spectrum at a range of redshifts, we can immedi- 
ately find angular diameter distance c/a and Hubble parameter 
H as a function of redshift. Having good knowledge of these 
dependencies allows us to put constraints on the properties of 
the dark energy. To carry out this project one needs a tracer 
population of objects whose clustering properties with respect 
to the underlying matter distribution is reasonably well un- 
derstood. There have been several works discussing the usage 
of galaxies (Blake & Glazebrook 2003; Hu & Haiman 2003; 
Linder 2003; Seo & Eisenstein 2003) and clusters of galaxies 
(Hu & Haiman 2003; Majumdar & Mohr 2004; Hiitsi 2005) for 
this purpose. What is most important is that already cuiTently 



existing galaxy redshift surveys have lead to the detection of 
acoustic features in the spatial distribution of galaxies, this way 
providing clearest support for the standard gravitational insta- 
bility picture of the cosmic structure formation. In the paper by 
Eisenstein et al. (2005) the detection of the acoustic "bump" in 
the two-point redshift-space correlation function of the SDSS 
' LRG sample is announced. The discovery of similar fea- 
tures in the power spectrum of 2dF ^ galaxies is presented in 
Cole et al. (2005). These results clearly demonstrate the great 
promise of the future dedicated galaxy redshift surveys like 
K.A.O.S.'' Similarly, useful measurements of the acoustic scale 
can be hoped by the planned SZ cluster surveys like the ones 
carried out by the PLANCK Surveyor * spacecraft and SPT ^ 
(Hiitsi 2005) and also with a large future photometric redshift 
surveys (Blake & Bridle 2005). In this paper we calculate the 
redshift-space power spectrum of the SDSS LRG sample find- 
ing evidence for the acoustic oscillations down to the scales of 
~ 0.25 h Mpc"\ which effectively correspond up to the 9. peak 
in the CMB angular power spectrum. These scales in the CMB 
are very strongly damped due to the finite width of the last- 
scattering surface and also due to the Silk damping (Silk 1968). 
Also, at those scales the secondary CMB anisotropies (mostly 
thermal Sunyaev-Zeldovich effect (Sunyaev & Zeldovich 1972, 
1980)) start to dominate over the primary signal. On the other 
hand, features in the matter power spectrum, although being 
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Fig. 1. Comoving number density of galaxies as a function of comov- 
ing distance. Smooth solid line shows a cubic spline fit to the number 
density estimated for 50 discrete radial bins. 



small (~ 5% fluctuations), are preseved by the linear evolu- 
tion and so opening up the way to probe acoustic phenomena 
at scales smaller than the ones accessible for the CMB studies. 

The paper is structured as follows. In Sec. 2 we describe the 
dataset to be analysed, Sec. 3 presents the method of the power 
spectrum calculation. The main results of this work are given 
in Sec. 4 and finally we conclude with Sec. 5. 



2. Data 

We analyze the publicly available data from the SDSS Data 
Release 3 (Abazajian et al. 2005). Specifically, we carry out our 
power spectrum measurements using the subset of the SDSS 
spectroscopic sample known as the Luminous Red Galaxy 
(LRG) sample. The LRG selection algorithm (Eisenstein et al. 
2001) selects ~ 12 galaxies per square degree meeting specific 
colour and magnitude criteria The resulting set of galax- 
ies consists mostly of an early types populating dense cluster 
environments and as such are significantly biased (bias factor 
b ~ 2) with respect to the underlying matter distribution. The 
selection method is very eff'ective producing a galaxy sample 
with a reasonably high density up to the redshift of z ~ 0.5. 

Since the selection criteria are very complicated, involving 
both cuts in magnitude and in color, and also due to the steep- 
ness of the luminosity function the usual method of using only 
the luminosity function to determine radial selection function 
does not work here (Zehavi et al. 2005). Here we simply build 
the radial selection function as a smooth spline fit to the num- 
ber density profile as shown in Fig. 1 . To calculate distances we 
choose the cosmological parameters as given by the WMAP ^ 
'concordance' model (Spergel et al. 2003). Unfortunately we 
could not find the data describing properly the survey geom- 
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Fig. 2. Angular distribution of galaxies given in the SDSS survey co- 
ordinates (A, Tj). The survey mask is shown with solid lines. 



etry and so were forced to build the angular survey masks 
using the galaxy data itself. As the number density of galax- 
ies in the sample is rather high, one can determine relatively 
accurately the beginning, ending and also possible gaps in the 
scan stripes. The angular distribution of the galaxies and also 
the boundaries of the survey mask built in the above mentioned 
way is shown in Fig. 2. Here the angular positions are plotted 
using the so-called survey coordinate system of the SDSS ^. 
We apply lower and upper redshift cutoffs of 0.16 and 0.47 as 
also done in Eisenstein et al. (2005). In the analysis presented 
in this paper we have excluded the three southern stripes since 
these just make the survey window more anisotropic. We have 
also caiTied out calculations by including these stripes, finding 
out that the amount of additional information they provide is 
rather minor. Due to the very small sky coverage we have not 
used the data in slices 15, 42 and 43. So in total the analyzed 
galaxy sample covers 0.57 h^^ Gpc^ over 2920 square degrees 
on the sky and contains 37, 998 galaxies. 



3. Power spectrum calculation 

We calculate the power spectrum using a direct Fourier method 
as described in Feldman et al. (1994) (FKP). Strictly speak- 
ing, power spectra determined this way are the so-called pseudo 
spectra, meaning that the estimates derived are convolved with 
a survey window. Since in the case of the analysed LRG sample 
the volume covered is very large, reaching 0.57 h^^ Gpc^, and 
also the survey volume has relatively large dimensions along all 
perpendicular directions, the correlations in the Fourier space 
are rather compact. In general FKP estimator is working fine on 
intermediate scales and in the case the power spectrum binning 
is wide enough. More precisely, one can say that this estimator 
is fine for wavenumbers k and power spectrum bins larger 
than the width of the survey window function |WkP, where Wk 



* For the exact details of the selection criteria see Eisenstein et al. 
(2001) 



' http://lambda.gsfc.nasa.gov/product/map/ 



* The data related to the geometry seemed to be newer than the 
actual galaxy data. 

' The transformations between various coordinate systems used by 
the SDSS are given e.g. in Stoughton et al. (2002). 
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Fig. 3. Angle-averaged survey window |VVj.|-. The gray shaded area 
shows the region where the window is larger than 1% of its maximum 
value. 

is the Fourier transform of the product of the survey mask and 
radial selection function. The angle average of this quantity is 
shown in Fig. 3, where we have also marked with a gray shaded 
stripe the scale where it is larger than 1 % of its maximum value. 
This scale should serve as a rather conservative estimate of the 
effective window width. In the following we make sure that the 
power spectrum bins are chosen to be large enough so that the 
above discussed criteria are satisfied. To speed up the calcu- 
lations of the power spectrum we make use of the Fast Fourier 
Transforms (FFTs), meaning that we need to build density field 
on a grid first. We use 512^ grid with a Triangular Shaped 
Cloud (TSC) mass assignment scheme (Hockney & Eastwood 
1988). Very briefly the steps involved in calculating the power 
spectrum are as follows: 

1. determination of the survey selection function (including 
the survey geometry) and mean underlying number density 
n, 

2. calculation of the overdensity field on a grid using TSC 
mass assignment scheme, 

3. Fourier transformation of the gridded density field, 

4. calculation of the raw 3D power spectrum, 

5. subtraction of the shot noise component from the raw spec- 
trum, 

6. isotropization of the shot noise corrected 3D spectrum i.e. 
averaging over ^-space shells, 

7. application of the normalization correction due to selection 
effects, 

8. deconvolving the smearing effect of the TSC mass assign- 
ment. 

More detailed description of our computational method 
with the results of application to a multitude of test problems 
can be found in Hiitsi (2005). The power spectrum eiTors are 
estimated using the simple "mode counting" result of FKP (see 
also Tegmark et al. (1998)): 




k[hMpc''] 

Fig. 4. Power spectrum of the SDSS LRG sample with the bin width 
Ak ~ 0.005 h Mpc"' . The solid line shows the linearly evolved matter 
power spectrum for the "concordance" cosmological model multiplied 
by the square of the bias parameter b = 1.95. 

where Vk = Ank^Ak/iln)^ is the volume of the fc-space shell 
and Veff is the effective volume given by: 



Here n{z) is the mean number density of objects at redshift z, 
dVc is a comoving volume element and the weight function: 



w(z) 



const 
h(z) 

n(z) 
l+n(z)/' 



for volume weighting 
for number weighting 
for an optimal FKP weighting 



AP 
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In the following we use the FKP weight function, although due 
to the rather high number density and clustering strength of the 
sample the pure volume weighting is giving almost identical 
results. 



4. Results 

In this section we present the main results of this work. In Fig. 
4 we show the redshift-space power spectrum of the above 
described SDSS LRG sample. We see that the power spec- 
trum agrees remarkably well with the linearly evolved "con- 
cordance" model spectrum (shown with a solid line) down 
to the scales of ~ 0.1/zMpc beyond which the nonlin- 
ear evolution starts to boost the power. Here the bin width 
Ak ~ 0.005 h Mpc"' is probably too narrow for the bins to be 
independent of each other. In Fig. 5 we have increased the bin 
width up to ~ 0.02 /z Mpc"' and built a smooth cubic spline 
curve through the data points, revealing clear hints for the os- 
cillatory behaviour of the spectrum. In addition we have plotted 
the same linearly evolved theoretical power spectrum as in the 
previous figure as well as the spectrum of a similar model but 
with a zero baryon content (shown by the dashed curve with- 
out any oscillatory behaviour). In order to display the detected 
fluctuations in a more clear manner we have fitted a smooth 
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Fig. 5. Power spectrum of the SDSS LRG sample with the bin width 
AA: ~ 0.02 /iMpc"'. The upper solid line shows the smooth cubic 
spline model built through the data points, whereas the lower one cor- 
responds to the same theoretical linearly evolved model as in Fig. 4. 
The dashed line displays the model with zero baryon content. 

function in the form of a second order polynomial in log-log 
coordinates to the observed bandpowers. In Fig. 6 we plot the 
same spectra as in the previous figure after dividing by this 
smooth function. The light gray shaded stripe shows the effec- 
tive width of the survey window as described above. 

It would be interesting to study how these oscillations in 
the power spectrum translate to the peak in the two-point cor- 
relation function discovered by Eisenstein et al. (2005). For this 
purpose we use the cubic spline model described above and ex- 
tend it outside of the observed range by smoothly joining it to 
the power spectrum of the "concordance" model. The correla- 
tion function is now simply calculated as a Fourier transform of 
the power spectrum. To study the significance of the features at 
higher wavenumbers we have calculated correlation functions 
for several models that have oscillations "switched off" at var- 
ious scales. The spectra of these models are shown in Fig. 7, 
where for the shake of clarity we have introduced slight verti- 
cal shifts between the curves so that the scales where the tran- 
sition to the featureless spectrum takes place are easily visible. 
The coiTesponding correlation functions are given in Fig. 8. As 
expected, we see how the peak in the correlation function is get- 
ting broader and also decreasing in amplitude as we erase more 
features in the power spectrum. This clearly demonstrates the 
importance of many of the up-downs in the power spectrum to 
produce a relatively sharp feature in the two-point correlation 
function. Additionally, in this figure the circles with errorbars 
show the correlation function determined straight from the data 
using the edge-corrected estimator given by Landy & Szalay 
(1993): 



DD - 2DR + RR 
RR ' 



which has minimal variance for a Poisson process. Here DD, 
DR and RR represent the respective normalized data-data, data- 

With the only exception that we have not plotted a spectrum for a 
model with zero baryon content. 



random and random-random pair counts in a given distance 
range. Random catalogs were generated with 25 times the num- 
ber of objects in the main catalog. The displayed errorbars are 
the simple Poissonian ones multiplied by the factor of 2. Since 
our aim here is not to perform any detailed correlation function 
study, this approximate treatment of errors, where one simply 
multiplies the Poissonian errorbars (which underestimate true 
errors) by some "fudge" factor in the range 1.5 - 2.5 (see e.g. 
Mo et al. (1992); Martinez et al. (1993)), should be well justi- 
fied. The crosses with the dashed-line errorbars show the two- 
point function as determined by Eisenstein et al. (2005). We 
see that in general our results agree very well with their calcu- 
lations. 

To estimate the period of oscillations we calculate the peri- 
odogram Pxioj) as described in Horne & Baliunas (1986): 



Pxito) = - 



where 



tan(2wT) - 



No 

2 X(tj) COS to(tj - t) 



No 

2 X(t j) sin Lo(t J -t) 



No 

2 COS^ CL>(tj - t) 



No 

2 sm2(Ajtj 



No 

sin LL>(tj - t) 

7=1 



No 

2 cos2wfj 



and X{tj) represents a time series with A^o data points. We cal- 
culate the frequency of oscillations for two different cases: (i) 
points with wavenumbers up to ^ ~ 0.2/iMpc"' (first 10 data 
points), (ii) all 15 data points. The periodograms normalized 
by the total variance Piv((y) - fx(w)/cr^ are shown in Fig. 9, 
where the solid line corresponds to the case (i). The uncertainty 
in the oscillation frequency is estimated as (Horne & Baliunas 
1986): 



37rcr^ 



where cr^ is the variance of the residuals after the signal has 
been subtracted, T is the total length of the data set and A is the 
amplitude of the signal. 

The acoustic scales measured in the above described man- 
ner are (107.8 + A3)h-^U^c and (101.8 ± 3.2)h-^U^c for 
the cases (i) and (ii), respectively. These should be compared 
with the best-fit WMAP "concordance" model prediction of 
106.5 /i ' Mpc. It is worth pointing out that the acoustic scale 
in the matter power spectrum is slightly larger than the corre- 
sponding scale in the CMB angular power spectrum due to the 
fact that acoustic waves do not stall completely at the epoch 
of decoupling, but keep on moving (although with a strongly 
reduced speed) up to the redshift z ~ 1 00 (the so-called drag- 
epoch). 

5. Discussion and Conclusions 

In this paper we have calculated the redshift-space power spec- 
trum of the SDSS LRG sample, finding evidence for a series of 
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Fig. 6. The same spectra as in Fig. 5, except divided by the best fit- 
ting parabola in log-log coordinates, revealing clearly the oscillatory 
behaviour. The gray shaded stripe shows the efi'ective width of the 
survey window as described in the main text. 
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Fig. 7. Model spectra used to study the significance of the oscillatory 
features. Slight vertical shifts have been introduced in order to show 
clearly the locations of the transition scales applied. 

acoustic features down to the scales of ~ 0.25 /iMpc"'. Using 
the obtained power spectrum we predict the shape of the spa- 
tial two-point coiTelation function, which agrees very well with 
the one obtained directly from the data. Also, the directly cal- 
culated correlation function is consistent with the results ob- 
tained by Eisenstein et al. (2005). We have made no attempts 
to put constraints on the cosmological parameters, rather we 
have assumed in our analysis the "concordance" cosmologi- 
cal model. The derived acoustic scale (107.8 + 4.3)/!"'Mpc 
(k < 0.2/!Mpc"') agrees well with the best-fit WMAP "con- 
cordance" model prediction of 106.5/; 'Mpc. Using all the 
data down to the wavenumbers k ~ 0.3 /iMpc the obtained 
acoustic scale is somewhat shorter, which might be due to the 
nonlinear effects (transfer of power from large to small scales). 

The existence of the baryonic features in the galaxy power 
spectrum is very important, allowing one (in principle) to ob- 
tain Hubble parameter H and angular diameter distance as 
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Fig. 8. Two-point correlation functions corresponding to the models 
shown in Fig. 7. The circles display the correlation function as deter- 
mined directly from the data. The plotted errorbars are simply 2 times 
the Poissonian errors. And finally, the crosses with dashed errorbars 
show the correlation function found by Eisenstein et al. (2005). 



a function of redshift, this way opening up a possibility to con- 
strain properties of the dark energy (Hu & Haiman 2003). The 
currently existing biggest redshift surveys, which are still quite 
shallow, do not yet provide enough information to carry out 
this project fully. On the other hand, it is extremely encour- 
aging that even with the current generation of redshift surveys 
we are already able to see the traces of acoustic oscillations 
in the galaxy power spectrum, showing the great promise for 
the dedicated future surveys like K.A.O.S. We have seen that 
acoustic features seem to survive at mildly nonlinear scales 
(k > 0.1 /iMpc '), which is in agreement with the results of 
the recent N-body simulations (Springel et al. 2005; Seo & 
Eisenstein 2005). In order to fully exploit available information 
one needs a complete understanding of how nonlinear effects 
influence these features. Nonlinear bias and redshift space dis- 
tortions also add extra complications. In general redshift-space 
distortions, biasing and nonlinear evolution do not create any 
oscillatory modulation in the power spectrum and so acous- 
tic features should be readily observable. On the other hand, 
only the effects that change the amplitude of the spectrum (e.g. 
redshift-space distortions on large scales) are easily factored 
out, but the ones involving stretching/compressing of the spa- 
tial scales do start to interfere with the cosmological distor- 
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Fig. 9. The total variance normalized periodograms Pn{Q = 
Px(L)/cr^ for (i) power spectrum points with wavenumbers up to 
k ~ 0.2/iMpc"' (first 10 data points), (ii) all 15 data points. Case 
(i) is shown with a solid line. 



tions and so inhibit our ability to draw conclusions about the 
underlying cosmological model. So far there have been only a 
few works studying these important issues (e.g. Springel et al. 
(2005); Seo & Eisenstein (2005); White (2005)) and probably 
it is fair to say that currently we really do not have a full theo- 
retical description of them. This is also the reason why we have 
not tried to carry out any parameter estimation in this paper. 

The bare existence of the baryonic oscillations in the galaxy 
power spectrum tells us something important about the under- 
lying cosmological model and the mechanism of the structure 
formation. First, it confirms the generic picture of the gravi- 
tational instability theory where the structure in the Universe 
is believed to be formed by the gravitational amplification of 
the small perturbations in the early Universe. Under the linear 
gravitational evolution all the density fluctuation modes evolve 
independently i.e. all the features in the power spectrum will 
be preserved. And certainly, we are able to identify features in 
the low redshift galaxy power spectrum that correspond to the 
fluctuations seen in the CMB angular power spectrum (which 
probes redshifts z ~ 1100), providing strong support for the 
above described standard picture of the structure formation. 
Actually, we can also probe scales that are unaccessible for 
the CMB studies due to the strong damping effects and steeply 
rising influence of the secondary anisotropics, reaching effec- 
tively the wavenumbers that correspond to the 8.-9. peaks in 
the CMB angular power spectrum. Second, the ability to ob- 
serve baryonic features in the low redshift galaxy power spec- 
trum demands rather high baryonic to total matter density ratio. 
In Blanchard et al. (2003) it has been shown that it is possible to 
fit a large body of observational data with an Einstein-de Sitter 
type model if one adopts low value for the Hubble parameter 
and relaxes the usual assumptions about the single power law 
initial spectrum. In the light of the results obtained in our paper 
these models are certainly disfavored due to the fact that the 
high dark matter density completely damps the baryonic fea- 
tures. And finally, purely baryonic models are also ruled out 



since for them the expected acoustic scale would be roughly 
two times larger than observed here " . So the data seems to 
demand a weakly interacting nonrelativistic matter component 
and all the models that try to replace this dark matter compo- 
nent with something else e.g. modifying the laws of gravity 
might have hard time to fit these new observational constraints. 
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